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(54) Title: APPARATUS AND METHOD FOR COMPRESSING INFORMATION 
(57) Abstract 

An apparatus and method for signal, image, or video compression 
chat achieves high compression efficiency in a computadonally efficient 
manner and corresponding decoder apparatus and methods arc disclosed. 
This technique uses zerotree coding (126) of wavelet coefficients (100) 
in a much more efficient manner dian the previous techniques. The 
key is the dynamic generaoon of the list of coefficient indices (114- 
120) to be scanned, whereby the dynamically generated list (114-120) 
only contains coefficient indices for which a symbol must be encoded. 
This is a dramatic improvement over the prior art in which a stadc list 
of coefficient indices is used and each coefficient must be individually 
checked to see whether a) a symbol must be encoded, or b) it is 
completely predictable. Addidonally, using dynamic list geneiadon (1 14- 
120), the greater the compression of the signal, the less time it takes co 
perform the compression. Thus, using dynamic list generanon (114- 
120), Che computational burden is proportional to the size of the output 
compressed bit stream (134) instead of being proportional to the size of 
the input signal or image. 
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APPARATUS AND METHOD FOR COMPRESSING INFORMATION 
The invention relates to an apparatus and method for efficiently 
compressing information for storage or transmission and an apparatus and 
method for decoding the compressed information. 
5 BACKGROUND OF THE INVENTION 

Data compression systems are useful to represent information as 
accurately as possible with a minimxmi number of bits and thus minimize the 
amount of data which must be stored or transmitted in an information 
storage or transmission system. One of the primary means of doing this is to 
10 remove redundant information from the original data. In the Proceedings of 
the International Conference on Acoustics, Speech and Signal Processing, 
San Francisco, CA March 1992, voltmie IV, pages 657-660, I disclosed a 
signal compression system which applies a hierarchical subband 
decomposition, or wavelet transform, followed by the hierarchical successive 
15 approximation entropy-coded quantizer incorporating zerotrees. The 
representation of signal data using a multiresolution hierarchical subband 
representation was disclosed by Burt et al in IEEE Trans, on Commun. Vol 
Com-31, No 4. April 1983, page 533. A wavelet pjrramid, or a critically 
sampled quadrature-mirror filter (QMF) subband representation is a specific 
20 kind of multiresolution hierarchical subband representation. A wavelet 
P3rramid was disclosed by Pentland et al in Proc. Data Compression 
Conference April 8-11, 1991, Snowbird, Utah. A QMF subband pjrramid has 
been described in "Subband Image Coding", J. W. Woods ed., Kluwer Academic 
Publishers, 1991 and I. Daubechies, Ten Lectures on Wavelets, Society for 
25 Industrial and Applied Mathematics (SIAM): Philadelphia, PA, 1992. 

Wavelet transforms otherwise known as hierarchical subband 
decomposition have recently been used for low bit rate image compression 
because it leads to a hierarchical multi-scale representation of the source 
image. Wavelet transforms are applied to an important aspect of low bit rate 
30 image coding: the coding of the binary map indicating the locations of the non- 
zero values, otherwise known as the significance map. Using scalar 
quantization followed by entropy coding, in order to achieve very low bit rates, 
i.e., less than 1 bit/pel, the probability of the most likely symbol after 
quantization - the zero symbol - must be extremely high. Typically, a large 
35 fi-action of the bit budget must be spent on encoding the significance map. It 
follows that a significant improvement in encoding the significance map 
translates into a significant improvement in the compression of information 
preparatory to storage or transmission. 
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To accomplish this task, a new data structure called a zerotree is 
defined. A wavelet coefficient is said to be insignificant with respect to a given 
threshold T if the coefficient has a magnitude less than or equal to T. The 
zerotree is based on the hypothesis that if a wavelet coefficient at a coarse 
5 scale is insignificant with respect to a given threshold T, then all wavelet 
coefficients of the same orientation in the same spatial location at finer scales 
are likely to be insignificant with respect to T. Empirical evidence suggests 
that this h3T)othesis is often true. 

More specifically, in a hierarchical subband system, with the exception 

10 of the highest fi-equency subbands, every coefficient at a given scale can be 
related to a set of coefficients at the next finer scale of similar orientation. 
The coefficient at the coarsest scale will be called the parent node, and all 
coefficients corresponding to the same spatial or temporal location at the next 
finer scale of similar orientation will be caUed child nodes. For a given parent 

15 node, the set of all coefficients at all finer scales of similar orientation 
corresponding to the same location are called descendants. Similarly, for a 
given child node, the set of coefficients at all coarser scales of similar 
orientation corresponding to the same location are called ancestors. With the 
exception of the lowest fi-equency subband, all parent nodes have four child 

20 nodes. For the lowest frequency subband, the parent-child relationship is 
defined such that each parent node has three child nodes. 

A scanning of the coefficients is performed in such a way that no child 
node is scanned before any of its parent nodes. Given a threshold level to 
determine whether or not a coefficient is significant, a node is said to be a 

25 ZEROTREE ROOT if 1) the coefficient has an insignificant magnitude, 2) the 
node is not the descendant of a root, i.e. it is not completely predictable fi-om a 
coarser scale, and 3) all of its descendants are insignificant. A ZEROTREE 
ROOT is encoded with a special symbol indicating that the insignificance of 
the coefficients at finer scales is completely predictable. To efficiently encode 

30 the binary significance map, three symbols are entropy coded: ZEROTREES, 
ISOLATED ZEROS, and non-zeros. 

SUMMARY OF THE INVENTION 
The invention is an apparatus and method for encoding information 
with a high degree of compression and corresponding decoder apparatus and 

3 5 methods. TlJ^invention uses zerotree coding of wavelet coefficients in a much 
more efficient manner than previous techniques. The key to this invention is 
the dynamic generation of the Ust of coefficient indices to be scanned, whereby 
the dynamically generated Ust only contains coefficient indices for which a 
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sjnnbol must be encoded. This is a dramatic improvement over the prior art 
in which a static list of coefficient indices is used and each coefficient must be 
individually checked to see whether a) a symbol must be encoded, or b) it is 
completely predictable. 
5 Apparatus for encoding information comprises filter means for forming 

a wavelet transform of the image, means for forming a zerotree map ft'om the 
output of the filter means, means for setting an initial threshold for encoding 
the wavelet coefficients, means for encoding the significant coefficients on an 
initial dominant list fi-om the coarsest level of the transform and the children 
10 of those coefficients whose indices are appended to the dominant list as the 
coefficient of the parent is found to be significant, means for reducing the 
threshold, means for refining the estimate of the value of the significant 
coefficients to increase the accuracy of the coded coefficients, and means for 
cycling back to means for encoding the significant coefficients to scan the 
15 dominant list anew at the new, reduced threshold. 

The invention is a method for encoding information comprising the 
steps of forming a wavelet transform of the image, forming a zerotree map of 
the wavelet coefficients, setting an initial threshold for encoding the wavelet 
coefficients, encoding the significant coefficients on an initial dominant list 
20 ft-om the coarsest level of the transform and the children of those coefficients 
whose indices are appended to the dominant list as the coefficient of the 
parent is foimd to be significant, reducing the threshold, refining the estimate 
of the value of the significant coefficients to increase the accuracy of the 
coded coefficients, and cycling back to scan the dominant list anew at the new, 
2 5 reduced threshold. 

The invention also includes an apparatus and method to decode signals 
that have been encoded using the apparatus of the invention. 

BRIEF DESCRIPTION OF THE DRAWING 
Fig. 1 is a schematic illustration of parent-child dependencies of 
30 subbands in an image decomposed to three scales. 

In Fig. 2, the parent-child relationship for three generations of a 
subsampled image is illustrated. 

Fig. 3 is a block diagram of the image encoder of Fig. 1 
Fig. 4 is a block diagram of the filter apparatus used to perform a one- 
35 scale subband decomposition of an image. 

Fig. 5, 6, and 7 are flowcharts illustrating the coding method of the 
invention. 

Fig. 8 is a block diagram of the image decoder of the invention. 
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Figs. 9, 10 and 11 are flowcharts iUustrating the decoding method of the 
invention. 

DETAILED DESCRIPTION 
For a wavelet hierarchical subband decomposition, the image is 
5 decomposed using times two subsampling into high horizontal - high vertical 
(HH) high horizontal - low vertical (HL), low horizontal - high vertical (LH). 
and low horizontal - low vertical (HH), frequency subbands. The LL subband 
is then further subsampled times two to produce a set of HH. HL, LH and LL 
subbands. This is down recursively to produce an array such as that 
10 illustrated in Fig. 1 where three subsampUngs have been used. Preferably six 
subsampUngs are used in practice. The parent-child dependencies are shown 
with the arrows pointing from the subband of the parent nodes to the 
subbands of the child nodes. The lowest frequency subband is the top left 
LLl, and the highest frequency subband is at the bottom right HH3. In this 
1 5 example, all child nodes hiave one parent. 

In Fig 2 the parent-child relationship for three generations of a 
subsampled image is illustrated. A single parent node 82 has four child nodes 
84 corresponding to the same region in the image with times four 
subsampling. Each child node 84 has four corresponding next generation child 
20 nodes 86 with a further times four subsampUng. 

In the previously disclosed method a Dominant List contains all 
coefficients for a given pass that have not yet been found to be significant. In 
the present invention, for a given pass, the Dominant List will contam the 
coordinates of those coefficients for which one of the following two conditions 
2 5 apply 

1 The coefficient has not yet been found to be significant AND, a symbol 
will be generated on the CURRENT dominant pass, or 

2 The coefficient has been previously found to be significant BUT one of 
its descendants satisfied condition 1). 

30 What is excluded from the Dominant List in the present invention is the 

coordinates of those coefficients which have not yet been found to be 
significant, but have ancestors that are zerotree roots. In this situation, no 
symbol would be generated. However, under the prior method, this position 
would have to be checked to verify this. The method of the invention avoids 

35 wasting computation in checking the predictabiUty of coefficients that are 
part of zerotrees. and thus predictably insignificant, by dynamically 
generating the Dominant List during the dominant pass. At the begmmng of a 
dominant pass, the Dominant List contains the coordinates of all of the 
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coefficients in the coarsest scale (lowest frequency ) subband whether or not 
they have been previously found to be significant. As this list is scanned 
during a dominant pass, if a coefficient in the coarsest scale subband has not 
previously been found to be significant, then one of 4 sjnnbols is encoded: 1) 
5 ISOLATED ZERO, 2) ZEROTREE ROOT, 3) POSITIVE SIGNIFICANT, 4) 
NEGATIVE SIGNIFICANT. Specifically, a ZEROTREE ROOT symbol is 
used when all of the descendants of the coefficient under consideration are 
either predictably insignificant, or have previously been found to be 
significant. In other words, a ZEROTREE ROOT implies that for all of the 

10 descendants, no symbol will be generated on the current dominant pass. In 
the current invention whenever a coefficient on the dominant list is not coded 
as a ZEROTREE ROOT, the Dominant List for the current pass is appended 
to include the children of the current coefficient. Thus, the Dominant List for 
the current pass will contain the coordinates of all coefficients that are not 

15 descendants of zerotrees, and thus the computational cost and compression 
time will be significantly faster than that of the prior art. 

Since the lowest frequency subband has children in three different 
subbands, during the scanning of the subband labeled LL3 in Fig. 1, three 
separate sublists of child coefficients to be scanned are maintained, one each 

20 for subbands labeled LH3, HL3, and HH3. Afler subband LL3 is finished, the 
three new sublists are combined and appended onto the dominant list. From 
this point on all coefficients have children in the same subband. Therefore, 
whenever a coefficient is not encoded as a ZEROTREE ROOT, the coordinates 
of its 4 children are appended to the end of the Dominant List. 

25 In the disclosed means for determining the scanning order of the 

wavelet coefficients comprises the encoder and decoder each maintain two 
lists of coefficient coordinates in memory. At the beginning of a dominant 
pass, a dominant list contains the coordinates of those coefficients that have 
not yet been found to be significant. A subordinate list contains the 

30 coordinates of those coefficients that have been previously found to be 
significant with respect to previous, larger, thresholds. Each list contains the 
coordinates in the order in which the coefficients will be further processed for 
coding. At the beginning of the first stage, all coefficient coordinates are found 
in an original dominant list and the subordinate list is empty, since no 

35 coefficient has yet been established as being significant. As coefficients are 
determined to be significant or insignificant in the dominant pass, their entries 
are moved fi-om the original dominant list to either the subordinate list, or to a 
newly generated dominant list, respectively. The original dominant list 
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20 



25 



disappears after the coefficients associated ^th all of its entries have been 
evaluated, and the new dominant list and the subordinate hst are 
subsequently refined at progressively finer thresholds. 

In that process, a dominant list contains all coefficients for a g.ven 
pass that have not yet been found to be significant. In the current invention 
for a given pass, the dominant Ust will contain the coordinates of those 
coefficients which meet at least one of the following thr« condifon : 1^ The 
coefficient has not yet been found to be significant AND, a symbol w.U be 
United on the CURRENT dominant pass, or 2) The coefficent has been 
previously found to be significant BUT one of its descendants satisfied 
condition 1) or 3) The coefficient is in the coarsest scale Oowest frequency) 

"""""^at is excluded in the preferred method as opposed to the previous 
„,tV,„H is the coordinates of those coefficients which have not yet been found 
The sl^^firnt but have ancestors that are ZEROTREE ROOTS. In this 
luaXrno symbol would be generated. -der «.e^rev,ously 

disclosed method, this position would have to be checked to venft' tlu8_^ 

■ This preferred method avoids wasting computation m checfang the 
predictability of coefficients that are part of zerotrees. and thus predictably 
C^cant^by dynamically generating the dominant list during the domman 
pasflt the beginning of a dominant pass, the dominant hst contains the 
coordinates of alTof the coefficients in the coarsest scale (lowest f^equencyj 
subband whether or not they have been previously found to be ^'S-^^"'- 
his ust is scanned during a dominant pass, if a coeffic^nt » 
scale subband has not previously been found to be -^^1^'^°^:^°'^^ 
symbols is encoded; 1) ISOLATED ZERO, 2) ZEROTREE KOO■I^ 3) 
P^ITIVE SIGNIFICANT, or 4) NEGATIVE SIGNIFICANT. SpecificaUyj^a 
ZEROTREE ROOT symbol is used when all of the ^-cendan's of *e 
coefficient under consideration are either predictably "-^^T^^;^',;; 
previously been found to be significant. In other words -^^^^ 
LpUes that for aU of the descendants, no symbol will be generated on the 
imphes tn preferred method whenever a current 

current dominani pass, f rrr:^Tyr\rrx>j?T? PHOT the 

coefficient on the dominant list is not coded as a ZEROTREE ROOT 'he 
Zlen of the current coefficient are appended to the dominant hst for *e 
c^Tnt pass Thus, the dominant list for the current pass wiU contain the 

cZltes of all co;fficients that are not ^--^.'^ "^^X, 
thus the computational cost and compression tune will be significantly faster 

than that of the prior art. 
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In the preferred mode of operation for image coding, since the lowest 
frequency subband has children in three different subbands shown in Fig. 1, 
during the scanning of the subband labeled LL3, three separate sublists of 
child coefficients to be scanned are maintained, one each for subbands labeled 
5 LH3, HL3, and HH3. After subband LL3 is finished, the three new sublists 
are combined and appended onto the dominant list. From this point on all 
coefficients have children in the same subband. Therefore, whenever a 
coefficient is not encoded as a ZEROTREE ROOT, the coordinates of its 4 
children are appended to the end of the dominant list. 

10 Fig. 3 is a block diagram of the image encoder 28 shown in Fig.l of the 

aforesaid concurrently-filed application Serial No. 08/007,465 for an image 
input 27. Image encoder 28 includes filter means 100 for the subband 
decomposition to produce the wavelet transform of the image 27. The wavelet 
transform output 102 of the filter means 100 provides the input to means 104 

15 for computing the zerotree map 106. The zerotree map is used to determine if 
an insignificant coefficient is a ZEROTREE ROOT or has at least one 
significant descendant. The zerotree map Z for each coefficient is formed as 
follows: 1) The magnitude of each coefficient in the subband decomposition is 
copied into the array for the zerotree map. 2) All bits in the copy except the 

20 most significant "1" bit for each coefficient are set to zero. This value is 
referred to as the initial zerotree map value. The actual zerotree map value for 
a coefficient is defined as the bitwise OR operation between the initial zerotree 
map value of itself and the initial zerotree map value of all of its descendants. 
This can be computed recursively, always computing the actual zerotree map 

25 value for children before their parents, by computing the bitwise OR operation 
between its initial zerotree map value and the actual zerotree map values of 
its children, which would have already been computed. Using the bitwise OR 
operation, a bit position in Z is "1" when either itself or any of its descendants 
has their most significant "1" bit at that position. 

30 The zerotree map 106 provides the input to means 108 for computing 

the initial threshold to be used in the analysis of the wavelet coefficients. The 
threshold is preferably set to values which are powers of two. The initial 
threshold is then preferably chosen to be a value equal to the largest power of 
two less than the largest zerotree map value in the coarsest LL subband. 

35 The output 110 of the means 108 provides the input to means 112 for 

encoding the significant coefficients on the initial dominant list firom the LL 
subband at the coarsest level (shown in detail in Fig. 5) and their children 
whose indices are appended to the dominant list as the coefficient of the 



PCTAJS94/00451 

WO 94/17492 .f- 

8 

. u o ,prntree root (shown in detail in Fig. 6), The 
parent is found not to be a zerotree root ^ ^^^^elet 
, N.;«T, 109 the zerotree map 10b, tne nvunuci 

116 for reducing the threshold, P^^^'f ^ subordinate list to 

116 provides the -J^Xiffi^nitho^ in detaiUn Fi. ^ The 

increase the accuracy of the coded coe ^^^^^ ^.^^^ 

encoding s^sten. then cdes ac. ~ ^th^shold. The output ^ of 
10 scan the dominant hst an w at *e ne^ ^ ^ ^^^^^ 

means 120 and an output 124 of ^„ f„ compression 

„aerl26. "^-^tr/:" t^^TrcI^ - Bata Comp^^^^^^^^^ 
is discussed m W.t en 3 ^ j^, iggv. and in Langdon. 

15 nTCarii^^^'— cCoding-.IBM,..^^ 

-^t'nrSrencodingthehitstreamh^er^^^^^^^^^ 
„^..r of wavelet sca^s 10. U^^^^^^^^^ 
.0 :dr2:"Zt of clr 12^ is the encoded, compressed hit stream 

representative of the input image_ decomposition of an 

I. Fig. 4 filter means 100 for a pass 
i^ge is shown. Input ^^^e \Xly before being horizontally 

filtered by filters 204 and /'^T ,„bjampling filters 208. The 
25 subsampled by times two h""^""^' vertically low 

horizontally low pass output ^^^f^'t^/^^rbefore being vertically 

pass and high pass filtered by filters a"d 21 ,,y. 

Lbsampled by times - --fir/p ^r^ortputs 216 and 218. 
Subband components LL and i.n pp ^^^^ ^^^^^ 

30 respectively. Similarly, the hon--'^^ ^^^^P-^^*^^ ^rs 210 and 212 
206 and 208 is vertically low subsampUng filters 

before being vertically subsampled "J '""^'^ outputs 220 

214. respective.. ^''''-VZTrorteTtith t ^apparatus 200 are 
and 222. respectively. Tlie fil er ^ horizontal 

35 preferably digital quadrature mi^or bands, 
and vertical fi-equency bands into low f J'"^^/^ „„,.,^le subband 

Additional scales are P"^-^ Q^t -H^^ 
decomposition on the LL subband. QMF mters 
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are similar to one another. A reference for this is "Subband Image Coding", 
edited by John W. Woods, Kluwer Academic Publishers, Boston, 1991 and I. 
Daubechies, Ten Lectures on Wavelets, Society for Industrial and Applied 
Mathematics (SIAM): Philadelphia, PA, 1992. More specifically, 
5 specifications of the low and high pass filters can be found in the Chapter 4, 
"Subband Transforms" of that book by E. P. Simoncelli et al. The Appendix of 
that Chapter provides details of the filter design. The decomposition itself is 
aliso disclosed in "A Practical Approach to Fractal-Based Image Compression" 
by A. Pentland and B. Horowitz, Proceeding of the Data Compression 

10 Conference 1991, Snowbird UT, page 179, Fig. 2. 

The bit budget may be a function of a variety of factors, such as the 
product of the desired bit rate (e.g., 0.25 bits/pel) and the nxmiber of 
coefficients to be coded (in this case 262,144 subband coefficients for a 512 x 
512 image), the capacity of a modem by which the coded bit stream will be 

15 transmitted, or the capacity of a disk or tape used to store the coded bit 
stream. 

The threshold is typically progressively reduced by a factor of 2, 
although other factors may be selected depending on the reqxiirements of a 
particular system. The threshold has to be a power of two for the bitwise 

20 operations to work and is preferred. In the dominant pass three coefficient 
possibilities (positive significant, negative significant and insignificant) are 
coded using an alphabet of four symbols. The four symbols are "0" 
representing an insignificant isolated zero coefficient, "1" representing an 
insignificant ZEROTREE ROOT, '*2" representing a positive significant 

25 coefficient, and "3" representing a negative significant coefficient. The '*1" 
ZEROTREE ROOT symbol is the only "grouped" symbol and can occur only in 
the case of a coefficient with descendants. Although five coding possibilities 
exist, only the four listed coding symbols are used. The uncoded fifth 
possibility, an insignificant coefficient descending fi-om a ZEROTREE ROOT, 

30 is implicitly coded when a "1" ZEROTREE ROOT sjrmbol is encoded for the 
ancestor. 

The method for for determining the significance of the coefficients fi-om 
the LL subband at the coarsest frequency on the initial dominant list and of 
the descen<^|^ts in all other subbands after they are appended to the 
35 dominant list are described with reference to Figs. 5 and 6, respectively. 
These operations are repeated in a loop for later, smaller thresholds. Coding 
starts by obtaining the coordinates of the first coefficients fi'om the low-low 
dominant list. The list is defined such that each entry is an ordered pair (x,y) 
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coefficient ha b en^^^^^^ ^-"^'V' 

T ^th a sTcaUed Peano-Hilbert plane fUUng curve as discussed m 
accordance -'^ aj caUed p^^^.. TransacUons 

' ^"^:rL . --r. «8e. . ri. U » dia^a^ or an 

.,.en..rs<^n.p^^^^^^^^^^ 

In F.g. 5 ^ Significance Flag bit 708. At system 

sign bit 704. ^^^'I'^'Zt^Z^ Flag bit 708 for all coefficients 
10 startup pnor U. ""^^^^^ ^dent is found to be significant using 
is initially set to FALSE (0). ^len a ^it will be set to 

,Ue means described in Figs. 5 and ; ^ J^, ^, 3„3ller 
TRUE (1) and wUl remain TRUE (1) lor all lu 

thresholds. dominant list includes 

15 In Fig. 5 at the start of c«^°8 J (coarsest) 

only the coo~-f ^^f : r^tf s^:^ with the step 710 of fetching 
subband eg LL3 n Flg_ 1. I^e - ^^^.^^^ 

the coordinates of the first coemc dominant list. 

712 of fetching the value x °f ^^^^ ttaeshold T. If the 

20 Step 714 Si^rra^t m« bit for the coefficient is checked in 

answer is NO, then the Significance nag ,„„tree map value Z is 

step 716. If the Significance Lue T in a 

fetched in step 71 and compared m su^^^^ ^^^^^^^ ^ ^ 

bitwise AND ^^•'^''^''^J^^Zb.O. If it is a ZEROTREE ROOT. 
25 ZEROTREE BOOT or an ISOLATED itB ^^^^ 

then a ™'^««°l;^^:lrrdit hel updated in step 724 and 
arithmetic coder. '^^ ^"'^ , l.Zr-!28 730 and 732 to update the bit 
the system recycled « l^''Z,..r.^ list, compare indices to 
budget, increment dominant list has been fiUly scanned and 

30 determine whether the If it is a ISOLATED 

determine whether the bit budget has been exce d^^ I 

TGOT ATFD ZERO symbol is generatea in siep 
ZERO, then coder model is then updated in step 

in the arithmetic coder. Th ^thm ^^^^^ 
736, the indices of the thre":hlWr ^ ^ ^^^^^ „ 

35 coiresponding dominant Ust^n step 738^«^ J^^^ 

rjrrcr,r:;et^rsofthetbreechildrenofthecoefficientare 
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still appended to the corresponding dominant list in step 738 and the system 
recycles are before. 

If the answer is YES at step 714, the sign of the coefficient is 
determined in step 740 to be POSITIVE or NEGATIVE, In step 742 (shown 
5 in dashed box), depending upon whether the coefficient is POSITIVE or 
NEGATIVE, a POSITIVE or NEGATIVE symbol is generated, coded into the 
arithmetic coder, and the arithmetic coder model is updated for a POSITIVE 
or NEGATIVE symbol respectively. The magnitude of the coefficient x is 
then is then appended to the Subordinate List in step 744. In step 746 the 

10 Significance Flag is then set to be TRUE and the magnitude of the coefficient 
X is set equal to zero. The three children of the coefficient are then appended 
to the corresponding Dominant List in step 738. 

This procedure is repeated for each coefficient on the LL Dominant List 
imtil all coefficients are coded and the end of the LL Dominant List is detected 

15 in step 730. At that point the three lists of children of the LL Dominant List 
coefficient are then concatenated to form a single list appended to the 
Donainant List and the coding process continued as indicated by step 752 and 
as illustrated in Fig. 6, 

In Fig. 6 at the start of the coding process the dominant list includes 

20 the coordinates of the coefficients which are children of coefficients of the 
lowest frequency (coarsest) subband which were not foimd to be zerotree roots 
for the current threshold T. The method starts with the step 810 of fetching 
the coordinates of the first coefficient on the child dominant list and the step 
812 of fetching the value x of the first coefficient on the child dominant list, 

25 Step 814 is comparing the magnitude x with the initial threshold T. If the 
answer is NO, i.e. the magnitude of x is less than T, then the Significance Flag 
bit for the coefficient is checked in step 816. If the flag bit is found to be 
TRUE, the coefficient is known to have been fotmd to be significant during a 
prior donainant pass and its children are appended to the Dominant List. If 

30 the flag bit is found to be FALSE then the zerotree map value Z is fetched in 
step 818 and compared in step 820 with the threshold value T in a bitwise 
AND (&) operation to determine whether the coefficient is a ZEROTREE 
ROOT or an ISOLATED ZERO. In step 822 (shown in dashed box), depending 
upon whether the coefficient is a ZEROTREE ROOT or an ISOLATED ZERO, 

35 a ZEROTREE ROOT or an ISOLATED ZERO symbol is generated, coded into 
the arithmetic coder, and the arithmetic coder model is updated for a 
ZEROTREE ROOT or an ISOLATED ZERO symbol respectively. If the 
coefficient is found to be a ZEROTREE ROOT, then the system recycled 
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«9R 830 and 832 to update the bit budget, increment the 
throughsteps826 828 830 and P^^^^ ^ ^^^^^^^ ^^^^^^^ 

VuJZZm^^'^^ and detennine whether the bit budget 
Donunant Ust »^ been mi y ^ ^^^^ ^^^^ ^he indices of the 

. ^^;rnVrct:L"-^^ 

.='^P«3«-"^*^='^''";"tt"erri4, the sign of the coefficient is 
If the ~ ^J^"p;;i:^„WtIVE. In step 842 (sho^ 
determined in step 840 to be FUaii coefficient is POSITIVE or 

in dashed box), depe^ng rJ^'^^toUst-^^^^ 
10 NEGATIVE, a U updated for a POSITIVE 

arithmetic coder, and the arithmetic ^ „f the coefficient is 

cr NEGATIVE symbol respectively, ^he -na^ude^o 

then appended to ^ t''— 

. ^'!:'r:u:r r.:: if- coefficient are then appended to the 

corresponding Dominant List in step 838 _ comlnant List 

This procedure is repeaud f r^^^^^^^^ 

TalTt^feTd of tt^oi^nalt Ust is detected in step 830. At that 
are coded and the end oi tne threshold level T the process 

.0 PO-^'-coeffi^entsba^^^^^^^^^^^ 

;rr:^:"^-tep8S0,.dasin..a^^^^^^^ 

the precision of the coefficients have tnus -^tiirns to the LL 

once the Subordinate Ust^^ comp e^^^^^^^ 

scanning as indicatedby sup 85^ ^.^^ ,U appended 

coefficlrh:rbe^col;rd^heSubordinateListcoeffi^^^^^ 

subordinate List during the scanofU^Bo^t^^^^^ 
is carried oat at a threshold level i„ , Mnary 

the threshold level is clearing the Botto. 

35 code by a st?fiig of ones and zeroes ine f 
Ust,tobeaefinedlat^A^^Mjs^^^^^^^^^ 
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906 initializes a reference bit string W equal to the bitwise AND of the 
magnitude x of the first coefficient on the Subordinate List and the mask M. 
Step 908 is the first step in the loop to scan all the coefficients on the 
Subordinate List for all previous thresholds. Step 910 sets a second reference 
5 bit string Y equal to the bitwise AND of the magnitude x of the first coefficient 
on the Subordinate List and the mask M. This step is repetitive of step 906 
for the first coefficient but takes on different values as other coefficients are 
scanned. Step 912 compares the values of the first and second reference bit 
strings W and Y, The answer is YES for all coefficient on the Subordinate List 

10 on the first pass since all are greater that the first threshold level. 

Step 914 is a bitwise AND comparison of the magnitude x of the 
coefficient with the present threshold. For example, if in the binary 
representation x = 1????? and the threshold T = 010000 the YES/NO decision 
depends upon whether the second bit in x is a one or zero. If a one, then the 

15 magnitude of x is greater than a level of 11???? or in the upper half of the range 
between threshold T and the previous threshold 2T. If a zero, then the 
magnitude of x is less than a level of 11???? or in the lower half of the range 
between threshold T and the previous threshold 2T. Thus an additional bit of 
accuracy has been gained in assessing the value of the particular coefficient. 

20 If a YES, the symbol "1" is coded into the arithmetic coder and the model is 
updated in steps 916 and 918, respectively. If a NO, the symbol "0" is coded 
into the arithmetic coder and the model is updated in steps 920 and 922, 
respectively. The coefficient is then removed firom the Subordinate List and 
added to a separate list, the Bottom List, which holds all coefficients in the 

25 lower half of the range. This is the sequence on the first scan of the 
Subordinate List since all coefficients are greater than T. 

The system then cycles through steps 926, 928, 930 and 932 which 
update the bit budget, increment the index on the Subordinate List, determine 
whether the end of the Subordinate List has been reached or not and tests to 

30 determine whether the bit budget for the storage or transmission channel has 
been exceeded. The system is then ready to test the next coefficient. 

Once all coefficients on the Subordinate List for that subordinate pass 
have been evaluated, the end of the Subordinate List has been reached and 
the system exits the scanning loop. Step 934 appends all coefficients on the 

35 Bottom List to the end of the Subordinate List. This step restores the 
coefficients on the Bottom List to the Subordinate List and places them at 
the end of the list, thereby completing a sorting the coefficients by their value. 
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Once the next supplemented by all those 

completed, the Subordxnate L.st has been 

5 coefficients having values between T ^-^^^ J^^ ^..^ M equal 
except when the bitwise AND of ^ fi„t .efficient on 

to Y does not equal the "AND of the mask M ^ 
the Subordinate list equal to W.certa^ to be a 1^ 

the entire .'''^'^'^^" tto^^List immediately in front of 

16 subordinate List is " * ,f:^;TB„ttom Ust and step 942 sets Y 

rw^"^tsr:^sfr::istothe.^^^^^^^^^^ 
::rr:m!:i:rrrre:i= 

15 factor of two from the previous threshold 

The decoder 1000, which operates m ttie conve 
Fi.. 3. includes means 1002 for ^^^f^^^^^Zuu^r^^Zi^ as an 
separate the header bits from '"^^ f^'^'"'^^^^ ^^Zn^ Juts to decoder 
input 1006 to the Dominant List f ^f ^j^^^^-^^^'^teshold 1025 and the 

,0 1008 are the image l^ouU of the means 1002. 

number of wavelet scales 1026 "bicn " 

The operation of decoder 1008 "'^^^jjfoos aTthe Dominant Ust 
1010 indicates the loop performed <^«"^^^ dominant List, 

with the appended children are ^^^^^.^..m. the threshold is 

scan the Dominant List at the new ttoeshold. ^^^^ 
The outputs of the ^^t.^ co~ to means 1022 for 

respectively, provide '^"^^"^^"^''^.aditional inputs to means 1022 
constructing the wavelet transform ^y^ ^ ^^^^ the output 

are the image size 1024. the number "l^'^"^!^ ^^^^^.^ Rags to be 
1028 of means 1030 for initializing the by f ^'f „„^„t 1032 of 
35 FALSE and the magnitudes of all TepresentaLe of the 

means 1022 is the array of wavelet '^ff'"^'' ^ ^ 1034 for 

transmitted image. The output 1032 provides » -P^' 
computing the inverse wavelet transform of the array 
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Additional inputs to means 1034 are the image size 1024 and the number of 
wavelet scales 1026. The method for performing the inverse wavelet 
transform has been described in "Subband Image Coding", J, W. Woods ed., 
Kluwer Academic Publishers, 1991 and I. Daubechies, Ten Lectures on 
5 Wavelets, Society for Industrial and Applied Mathematics (SLAM): 
Philadelphia, PA, 1992. 

The output 1036 of the means 1034 provides an input to means 1038 
for adding the image mean to the inverse wavelet transform. The image mean 
output 1040 from decoder means 1002 provides a second input to means 
10 1038. The output 1042 means 1038 is the transmitted image reconstructed 
to the extent possible given the number of bits transmitted. 

Figs. 9-11 are flowcharts illustrating the decoding procedure for a signal 
compressed using the procedure described in Figs. 5-7. 

In Fig. 8 at the start of the decoding process the dominant list includes 
15 only the coordinates of the coefficients from the lowest frequency (coarsest) 
subband, eg LL3 in Fig. 1. The decoder knows the protocol for the order in 
which the coefficients are scanned on compression, so that it knows the 
location of all coefficients. The method starts with the step 1102 of fetching 
the coordinates of the first coefficient on the LL dominant list followed by the 
20 step 1104 of checking the identity of the Significance Flag. If the Flag is 
TRUE, ie the magnitude of the coefficient is greater than an initial threshold 
Tq, then the indices of the three children of the coefficient are appended to the 

dominant list in step 1106. The system is then cycled through steps 1122 and 
1124 to increment the index on the dominant list, compare indices to 

25 determine whether the initial LL dominant list has been fully scanned. In the 
first Dominant List pass, all Significance Flags have the initial value FALSE. 
However, once a Significance Flag is set to be TRUE, it will remain TRUE for 
the entire decoding process on all fiiture Dominant List passes. 

Step 1102 fetches the next coefficient in the scan. If the Significance 

30 Flag is not TRUE, i.e., the magnitude of the coefficient is less than 2T, then a 
symbol from the arithmetic decoder is decoded in step 1108. In step 1110 a 
comparison is made to determine if the end of the bit stream has been 
reached. If YES, the decoding process is complete. If NO, then steps 1112 
and 1114 in sequence determine whether the sign bit is POSITIVE or 

35 NEGATIVE. Step 1116 then equates the coefficient to either +3T/2 or -3T/2 
depending upon whether the sign bit is POSITIVE or NEGATIVE. The model 
is then updated for a POSITIVE or NEGATIVE symbol. The position of the 
coefficient is then appended to the Subordinate List in step 1118 and the 
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i^aioes of the three '^'^^^^^{ZT-^ done at the encoder in step 
corresponding dominant list m step uu increment 
„S. The system is then ..cied the 
5 the index on the dominant hst f°^P 

iniUal IX domin^t hat '■"^^ ^^j^""^ ^OATTVE. then the symhol is 
If the symbol - ^^j,, ^^^i f„, a ZEROTREE ROOT, 

tested in step 1126 to detemnne if it ^^01 symbol in step 

„ VES. then the model is ^^ f^^'X^^ZZel is updated for the 
10 1128 and the system "^^''^^ J „f ,he three children of 

ISOLATED ZERO symbol m Step U28^ .^ing dominant Ust in step 

the coeffident are each appended to the corr p 

1106, and the system the LL Dominant Ust 

This procedure is «P«''^* '^/^'■^ ^^^^^ ^ dominant List is detected 

;"^rL:::irep=r^:"^^^ 

step 1134 and as illustrated in Fig. 10. d„^„ant Ust includes 

,nFig.lOatthestartofthedecoto„ro-s^ 

the coordinates of the '^o'^''^"'^/"^"" ' ^.^ce to Fig. 9. The method 
during the decoding process ^^^^"'^l^'^;^;^^''Jl, Lt coefficient on 
.tarts with the step 1202 o«etching * ^r-^ate ^^^^^^^ 
the child Dominant List followed "yt^e sUp 1 ^^^^^^ 
25 the Significance Flag. Hthe .lag ^a previous dominant 

is greater than 2T and — /""^^^J^ „,fficient are each appended 

pass, then the indices of the four children of ttie c^ ^^^^^ ^ 

L the corresponding dominant Ust in step 1106 -w-^^^ ^^^^ 

step 738. The system is then ^^^^^ to determine 
30 increment the index on the — ^*j^'^'jX3^^ 

whether the initial LL dominant list „xt scan. If 

Step 1202 then fetches the ne^ ^^^"^^^^t is less than ^T. then 

the flag is not TRUE, i.e. the magnitude "^^^ ^^OS. In step 1210 a 
a symbol from the arithmetic decoder . d^ed -^^^^^ ^^^^^^ ^ 

35 comparison is made to determine If the en 1212 

.cached. If VES, the decoding process. s^m^ete J ^^^^^ 

and 1214 in sequence d-'--«/';^f;;,*;,"tt either .3T/2 or -3T/2 
NEGATIVE. Step 1216 then equates the coefBcent 
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depending upon whether the sign bit is POSITIVE or NEGATIVE. The model 
is then updated for a POSITIVE or NEGATIVE symbol. The position of the 
coefficient is then appended to the Subordinate List in step 1218 and the 
Significance Flag for this coefficient is then set to TRUE in step 1220. The 
5 indices qf the children of the coefficient are appended to the dominant list in 
step 1206. The system is then cycled through steps 1222 and 1224 to 
increment the index on the dominant list, compare indices to determine 
whether the Dominant List has been fully scanned and start the loop anew. 

If the symbol is neither POSITIVE or NEGATIVE, then the sjrmbol is 

10 tested in step 1226 to determine if it is the symbol for a ZEROTREE ROOT. 
If YES, then the model is updated for the ZEROTREE ROOT s3niibol in step 
1228 and the system recycled. If NO, then the model is updated for the 
ISOLATED ZERO symbol in step 1230. The indices of the children of the 
coefficient are appended to the dominant list in step 1206 and the system 

15 recycles.. 

This procedure is repeated for each coefficient on the child Dominant 
List and all that are appended dtiring the cycling of this process until all 
coefficients are decoded and the end of the Dominant List is detected in step 
1224. At that point all the coefficients have been decoded for that threshold 

20 level T. The process continues as indicated by step 1232 in which the 
threshold is cut in half and the system exits to the Subordinate List in step 
1234, as illustrated in Fig. 11, to increase the precision of the coefficients 
which have thus far been foimd to be significant. Once the scan Subordinate 
List is completed at that threshold, the system retmnis to the LL Dominant 

25 List and repeats the process illustrated in Figs. 9 and 10 at the new threshold 
T/2 until the end of the list is reached at the lowest threshold used. 

In Fig. 11, the refinement of the value of the coefficients placed on the 
Subordinate List have thus far been found to be significant is carried out at a 
threshold level of T/2. The first step 1302 is clearing the Bottom List, to be 

30 defined later. In step 1304 a reference value S is set equal to Y, the 
magnitude of the wavelet coefficient at the coordinates for the first element on 
the Subordinate List. In step 1306 the coordinates of the next coefficient on 
the Subordinate List are fetched. In step 1308, the magnitude Y for the 
coefficient is obtained from the Subordinate List. In step 1310, the value Y is 

35 compared to S. If equal, then the symbol from the arithmetic coder is decoded. 

In step 1312, the bit stream is then checked to determine if the end of 
the stream has been reached. If NO, the symbol is checked in step 1316 to 
determine if the symbol is a "1". If YES, then the model is updated for the "1" 
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3^.0. in step 1320 and T/2 is added - *e — ^^^^^^ - 
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otherwise. . , v-g descendants, there is a 

F„ each wavelet coeffi^en '^^^ Vourse of encoding dominant 
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zerotree map value and the current threshold are bitwise ANDed, which 
effectively polls the bit position in the zerotree map corresponding to the 
current threshold. If the result of the bitwise AND operation is non-zero, an 
ISOLATED ZERO symbol is encoded indicating the presence of at least one 
5 descendant that will first be found significant during the current dominant 
pass. If the resxilt of the bitwise AND operation is zero, then a ZEROTREE 
ROOT symbol is encoded indicating that no additional symbols are required for 
any descendants at the current threshold. 

The zerotree map is generated prior to any encoding as illustrated by 

10 means 104 for computing the zerotree map 106 in Fig. 3. In the generation of 
the zerotree map, there is one element corresponding to each wavelet 
coefficient, including those with no descendants. Once the entire zerotree map 
is generated, the memory used to hold the zerotree map values for the 
coefficients with no descendants can be de-allocated. 

15 The following operations are performed independently for each element 

individually: 

1 The magnitude of each wavelet coefficient in the wavelet transform 
(subband decomposition) is copied into the array for the zerotree map. 

2 All bits in the zerotree map's copy except the most significant "1" bit 
20 for each wavelet coefficient are set to zero. This causes each zerotree map 

value to be the largest power-of-two smaller than the magnitude of the 
corresponding coefficient. This value is referred to as the initial zerotree map 
value. It can be computed by Right Shifting (shifting all bits to the right in a 
binary representation of the value) until the value is 1 and then Left Shifting 
25 back to the correct position,where new bits added dxiring the shifts are always 
"0". 

The actual zerotree map value Z is defined as the bitwise OR operation 
between the initial zerotree map value of itself and the initial zerotree map 
value of all of its descendants. To determine this efficiently, the actual 

30 zerotree map value for a child is always determined before that of its parent. 
The actual zerotree map value of the parent is set equal to the bitwise OR of 
the initial zerotree map value of the parent and the actual zerotree map value 
of its children, which would have already been determined. Thus, this 
computation^gins in the highest frequency subband whose coefficients have 

35 children, i.e. LH2, HL2, and HH2 in Fig. 1, and continues to lower frequency 
subbands tmtil the entire zerotree map is complete. Thus, using the bitwise 
OR operation, a bit position in Z is "1" when either its corresponding coefficient 
or any of its descendants has their most significant "1" bit at that position. 
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compression schemes may be used. 



wo 94/17492 



PCT/US94/00451 



21 

IN THE CLAIMS: 

1. In a digital compression encoder responsive to an array of a relatively 
large number of multibit coefiBcients defining image data at a relatively high 
resolution that has been hierarchical subband decomposed by a wavelet 
5 transfoim at a given plurality of lower resolutions, whereby the lowest low 
(LL) subband comprises an array having the smallest nimiber of multibit 
coefficients of all the wavelet-transform subband arrays, each LL subband 
multibit coefficient comprises an ancestor of all multibit coefficients of all 
higher subbands descended from that LL subband multibit coefficient, and 

10 each of said multibit coefficients of said LL subband and said multibit 
coefficients of said descendents thereof includes a most significant bit; the 
improvement wherein said encoder comprises: 

means responsive to each multibit coefficient of said LL subband and 
said multibit coefficients of all of its descendents of said wavelet transform for 

15 deriving a zerotree map, wherein said zerotree map includes a zerotree-map 
multibit coefficient corresponding to each coefficient of said LL subband array, 
and each zerotree-map coefficient has a binary "1" value in a bit position 
thereof only if either its corresponding LL-subband-array coefficient itself or 
any of the descendents of its corresponding LL-subband-array coefficient has 

20 its most significant bit in that bit position; and 

coding means for deriving a data-compressed code, wherein said coding 
means starts with a list of coordinates of said relatively small number of 
coefficients of said LL subband of said array and a relatively high initial 
threshold value to sequentially derive said data-compressed code from the 

25 multibit coefficients of the wavelet- transform subband and the multibit 
coefficients of the zerotree-map first at said at said relatively high initial 
threshold value and then at one or more successively lower threshold values. 
2. The encoder of Claim 1, wherein: 

at least certain ordinally-arrranged bits of each multibit coefficient of 

30 said LL subband and each multibit coefficient of its descendents defines the 
magnitude of that coefficient, a first of said ordinally-arrranged bits 
corresponds to the magnitude of said most significant bit, and each other of of 
said ordinally-arrranged bits corresponds to one-half the magnitude of its 
immediately preceding bit; 

35 said irdtial threshold value is equal to one-half the largest zerotree-map 

value as represented by the most significant bit of all the bits of all the 
zerotree-map multibit coefficients; and 
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each of said successively lower threshold values is e,ual to one-half of 

its — ^-^'>^'=:^^ coding -neans con-pHses: 

3. The encoder otCtom !• successively responsive to each of 

during each cyde of operation thereof, ^^.j^ely responsive to each 

.econd cycUcally-operated lope means / J ^ j, first 

of suhordinate lisu of coordinates of n>utab.t CO— 

,0 cyclicany-operatedlo^c^ean^^^^^^ 

for deriving a ' J^f ^^'.^^ved dominant list is for use during a 
cycle of operation thereof which deri ^ „eans; 

subsecuent cycle of operation of said fi^;'^^^^^,,^^, , „„„ accurate 

operation of said second lo^c means; and ^^^^ 

said list of coordinates of said relauvely j„„i„ant list of 

20 said LL subband of said array "-""^"JJ^high initial threshold 

25 coefficient includes a sigmficance flag bit that .Sim 

and said first logic means compnsest ,g„i,,de Ix. of each 

means tor comparing, in turn the a ^^^es during the 

individual coefficient addressed hy^aid^nu^^"-^^ 
current cycle of operation of smdfctlopcm«u^ ^^^^_^ 

30 during that current ^^-^^rr ffide"'.^ s^crdinate Ust and then set 
append the magnitude X of that coefficient ^•'^^ sigmficance flag bit to 
t^^magnitude of that coefficient its'" employing the 
TRUE, (M. if isl S T and the significan^ "^.^'illcative ;f whether that 
--'^''VZ-^.tn:11z.t^ ^OOr. Whereby that 
- It:: :rS^^^ — codi. means to - .nerate . 
=^gT!°f^^nrtr : -O and generate a 
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ZEROTREE ROOT symbol code in response to said first logic means 
determining a coefficient to be a ZEROTREE ROOT, and (4) using the 
ISOLATED ZERO symbol code to update said model for an ISOLATED ZERO 
and using the ZEROTREE ROOT symbol code to update said model for a 
5 ZEROTREE ROOT; 

index appending means operative, unless said first logic means has 
determined that coefficient to be a ZEROTREE ROOT, following the 
comparison by said first logic means during the first cycle operation thereof of 
each individual coefficient for appending the index of each of the three children 

10 of that coefficient of the starting dominant list to a separate one of three 
corresponding new dominant lists, and following the comparison by said first 
logic means during the each subsequent cycle operation thereof of each 
individual coefficient for appending the index of the children of that coefficient 
of the current dominant list to a new dominant list for use during the next 

15 cycle of operation of said first logic means; and 

means effective after all the individual coefficients addressed by the 
starting dominant list of coordinates have been operated on during the first 
cycle of operation of said first logic means for concatenating the three 
corresponding new dominant lists into a single dominant list for use during the 

20 second cycle of operation of said first logic means, and effective after all the 
individual coefficients addressed by the current dominant list of coordinates 
have been operated on dxxring each subsequent cycle of operation of said first 
logic means for cutting the threshold value to the next lower one of said 
successively lower threshold values and initiating a cycle of operation of said 

25 second logic means. 

5. The encoder of Claim 3, wherein said second logic means comprises: 
means for establishing a bottom list of coordinates of multibit 

coefficients which is initially cleared; and 

initial means for determining the magnitude W by which a first 
30 coefficient on the subordinate list exceeds the current threshold value, 

means associated with each of said first and second logic means for 
updating said bit budget on each cycle of operation thereof and for terminating 
operation thereof if the bit budget is exceeded. 

6. In a digital compression decoder for reconstructing image data from a 
35 bit stream comprising data bits of arithmetic-coded symbols that define 

coefficient models for an array of multibit coefficients of wavelet-transformed 
hierarchical decomposed subbands of said image data; wherein said decoder 
has an a priori knowledge of the protocol for the order in which the wavelet- 
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transformed multibit coefficients of said array were scanned on compression 
so that a list of the coordinates of the multibit coefficients of the lowest low 
(LL) wavelet-transformed subband of said array are predetermined; the 
improvement wherein said decoder comprises: 
5 data decoding means responsive to said data bits applied thereto for 

deriving the data values of said wavelet-transformed multibit coefficients from 
said arithmetic-coded symbols that starts with a list of -ordinates o said 
relatively small number of multibit coefficients of said LL subband of said 
array and a relatively high initial threshold value to sequentially derive the 
10 valuL of said wavelet-transfonned multibit coefficients first at said relatively 
high initial threshold value and then at one or more successively lower 
threshold values. 

7 The decoder ofClaim 6, wherein: , ,r 
each of said successively lower threshold values is equal to one-half of 

15 its immediately preceding threshold value. 

8 The decoder of Claim 6, wherein said data decoding means comprises: 

■ setting means for initially setting the magnitude of each of the values of 
said wavelet-transformed multibit coefficients of said array to zero; 

arithmetic decoding means for decoding said arithmetic-coded symbols, 
20 first cycUcally-operated logic means successively responsive to each of 

dominant Usts of coordinates of multibit coefficients, certain decoded S3^bols 
and each of said threshold values for deriving a subordinate hst of coordinates 
of multibit coefficients and updating models of the values of jav^^^ 
transformed multibit coefficients of said array during each cycle of operation 

second cyclically-operated logic means successively resp^^e ^^ch 
of subrdinateLs of coordinates of multiHt coefficients f^^^ 
J, • ^= decoded sy^te'k, and a cccrrent one 

said threshoH v»lu«s f<. ^^^^ '^^^l „f wavelet-transformed 

30 coefficienta and "1>^=;°* "°^;'',^'"^t"cL of operation thereof during 
multibit coeffidents of sa^da^ayd^gja^^ P^^^^ 

each cyde of ope -anon ^''^'"^•"'"t.^y^y.™^ U^meaos; and 

value comprises the threshold value winch are usea awr^ 
operation of said first logic means. 
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9. The decoder of Claim 8, wherein each wavelet transform multibit 
coefScient of said array includes a significance flag bit that is initially set to 
signify FALSE by said setting means, and wherein said first logic means 
comprises: 

5 means responsive to the setting of said significance flag bit being 

FALSE for determining, in turn, for the decoded symbol corresponding to each 
individual coefficient addressed by said dominant list of coordinates during the 
current cycle of operation of said first logic means is an ISOLATED ZERO 
symbol, a ZEROTREE ROOT symbol or is neither; if neither, put a value 

10 equal to 3/2 of the threshold value T during that cycle in the corresponding 
^ coefficient of an output map, append coordinates of that coefficient to said 
subordinate list and set significance flag to TRUE; if an ISOLATED ZERO 
symbol, update model for ISOLATED ZERO symbol; and if a ZEROTREE 
ROOT symbol, update model for ZEROTREE ROOT symbol; 

15 index appending means operative only if said significance flag is TRUE 

or the decoded s3ncQbol has been determined to be other than a ZEROTREE 
ROOT symbol for appending the index of each of the three children of that 
coefficient of the starting dominant list to a separate one of three 
corresponding new dominant Usts, and following the comparison by said first 

20 logic means dxiring the each subsequent cycle operation thereof of each 
individual coefficient for appending the index of the children of that coefficient 
of the current dominant list to a new dominant list for use during the next 
cycle of operation of said first logic means; and 

means effective after all the individual coefficients addressed by the 

25 starting dominant list of coordinates have been operated on during the first 
cycle of operation of said first logic means for concatenating the three 
corresponding new dominant Usts into a single dominant list for use during the 
second cycle of operation of said first logic means, and eff^ective after all the 
individual coefficients addressed by the current dominant list of coordinates 

30 have been operated on during each subsequent cycle of operation of said first 
logic means for cutting the threshold value to the next lower one of said 
successively lower threshold values and initiating a cycle of operation of said 
second logic means. 

10. The decoder of Claim 6, wherein said bit stream comprises both said 
35 data bits and header bits defining image dimensions, said initial threshold 

value and the number of wavelet scales, and said decoding means further 
comprises: 
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data decoding means. ,.„tree map for use in digitally compres'-"" 

11 Ametbodforderivrnga^erotreemP 

encoding an a.ay of a -l^^]^^^ ^^I^^^t has .een hierarchical subhan. 
image data at a relatively bagh resd^ ^solutions 
decomposed by a ^.^^^d Jn^ri Jan array having the smallest 

whereby the lowest low OW subband ^^^^^.^^.f^^ ^ubband arrays 
number of multibit -'ffi"'"*^;' j^'^rfses an ancestor of all mvdub t 
each LL subband ~"'J^ed from that subband mulUb.t 

coefficients of all higher -^^-^*'^;^ents of said LL subband and -d 
coefficient. ^dea^^fs^;»'^^,^„,„,^ 
multibit coefficients of saia oe 

said method comprising ^'^P „^ttbit coefficient corresponding 

generating a separate 3uch that each zerotree-map 

eac^ coefficient of said LL ^^^^^^^^J^^,^ ,,ereof only if either >ts 
, coefficient has a binary "l" value » f'^^ or any of the descendents of 
::„.pondingI..subban^^^ 
its corresponding LL sudd 

that bit position. f„ generating a diptal c°»P«= 

12 A system comprising (1) sn e ^^^ber of multibit 

,5 code in response to an array resolution that has been 

coefficients defining image data at a J ^ ^i^,,^ 

hierarchical subband decomposed by « comprises an array 

^Cr resolutions, whereby the o»^^aX^_^^^_^^^ „^ ^„ „.ve.et- 

having the smallest number f """"^ ^^uMt coefficient comprises 
,n triform subband arrays, each LL -^^^^ ^^^^d, descended from 
' an ^tstor of all multibit c«fficen«^ all^i^ ^^^^^ ,^ai,ien.s of 
t^a^L subband multibit c""'^^';^,, „f descendents thereo 
r^d LL subband ^^dT^f rdecoder for decoding the digital 

said encoder includes (3) means ' ^ „f descendents of 

„,saidLLsubband^ds^ 
said wavelet transform for den g 
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includes a zerotree-map multibit coefBcient corresponding to each coefficient 
of said LL subband array, and each zerotree-map coefficient has a binary "1" 
value in a bit position thereof only if either its corresponding LL-subband- 
array coefficient itself or any of the descendents of its corresponding LL- 
5 subband-array coefficient has its most significant bit in that bit position, and 
(4) coding means for deriving a data-compressed bit stream comprising data 
bits of arithmetic-coded S3rmbols, wherein said coding means starts with a Ust 
of coordinates of said relatively small number of coefficients of said LL 
subband of said array and a relatively high initial threshold value to 

10 sequentially derive said data-compressed code from the multibit coefficients of 
the wavelet-transform subband and the multibit coefficients of the zerotree- 
map first at said at said relatively high initial threshold value and then at one 
or more successively lower threshold values; and 

said decoder includes (5) data decoding means responsive to said data 

15 bits of said arithmetic-coded symbols of said data-compressed bit stream 
applied thereto for deriving the data values of said wavelet-transformed 
multibit coefficients from said arithmetic-coded symbols that starts with a list 
of a priori known coordinates of said relatively small number of multibit 
coefficients of said LL subband of said array and a relatively high initial 

20 threshold value to sequentially derive the values of said wavelet-transformed 
multibit coefficients first at said relatively high initial threshold value and then 
at one or more successively lower threshold values. 
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